April 24, 2013 2:8 WSPC/INSTRUCTION FILE hadron'spectra'paper 



o 

u 

< 

m 
(N 



Ph. 



^ 



S 



International Journal of Modern Physics A 
© World Scientific Publishing Company 



Hadron spectra in p+p collisions at RHIC and LHC energies 



P. K. Khandai 

Department of Physics, Banaras Hindu University, Varanasi, 221 005, India 
pkkhandai&gmail. com 



P. Sett 
Nuclear Physics Division, Bhabha Atomic Research Center, Mumbai, 400085, India 

D ■ p. Shukla 



Nuclear Physics Division, Bhabha Atomic Research Center, Mumbai, 400085, India 

pshuklaSbarc.gov. in 



Department of Physics, Banaras Hindu University, Varanasi, 221 005, India 



> ■ V. Singh 

(N 
(N . 

\^ Received Day Month Year 

• , Revised Day Month Year 

o 

rf\ • We present the systematic analysis of transverse momentum (pr) spectra of identified 

hadrons in p+p collisions at RHIC {y/s = 62.4 and 200 GeV) and at LHC energies {y/s 
= 0.9, 2.76 and 7.0 TeV) using phenomenological fit functions. We review various forms 
of Hagedorn and Tsallis distributions and show their equivalence. We use Tsallis distri- 
Ly( ■ bution which successfully describes the spectra in p+p collisions using two parameters, 

KN ' Tsallis temperature T which governs the soft bulk spectra and power n which determines 

the initial production in partonic collisions. We obtain these parameters for pions, kaons 
and protons as a function of center of mass energy (y/s). It is found that the parameter 
T has a weak but decreasing trend with increasing y/s. The parameter n decreases with 
increasing y/s which shows that production of hadrons at higher energies are increasingly 
dominated by point like qq scatterings. Another important observation is with increasing 
y^, the separation between the powers for protons and pions narrows down hinting that 
the baryons and mesons are governed by same production process as one moves to the 
highest LHC energy. 

Keywords: Tsallis distribution; Hadron spectra. 

PACS numbers: 13.85.Ni, 24.85. +p 

1. Introduction 

The heavy ion collisions at Relativistic Heavy Ion Collider (RHIC) and Large 
Hadron Collider (LHC) are performed to study strongly interacting matter at high 
energy density.'i' The measurements in p + p collisions are important to understand 
particle production mechanisrrP and are also used as baseline for heavy ion (HI) 
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collisions. The hadron spectra provide insight into particle production as well as 
interaction in the hadronic and quark gluon plasma (QGP) phases. The hadron 
spectra at high transverse momentum (pt) arise from fragmentation of high pT 
partons or jets from initial hard partonic collisions and are known to follow a power 
law distribution. The high pT hadrons are very important for QGP studies as they 
measure jet quenching^! effect in QGP. The low pt hadrons form the bulk of the 
spectra arising from multiple scatterings and follow exponential distribution depic- 
tive of particle distribution in a thermal system which may be more meaningful 
for heavy ion system. In addition, for heavy ion systems the hadron spectra at 
intermediate pT can arise from quark recombination.'^'^ 

We review various forms of Hagedorn and Tsallis distributions which are used to 
describe p+p collisions at RHIC and LHC energies and show their equivalence. The 
Tsallis distributiorP'21 can reproduce the full spectral shape of hadrons with just 
two parameters; temperature T and a non-extensive parameter q which is related 
to a measure of temperature fluctuations and the degree of non-equilibrium in the 
system."^' The Tsallis functional form is mathematically same as Hagedorn distri- 
butiorP with the parameter q related to power n of Hagedorn distribution which 
governs initial partonic collisions.'^ The Hagedorn function has been successfully 
used to describe the meson spectra in p-l-p collisions at RHIC energies with an 
additional parameter to describe heavy ion collisions as well.'^ 

In p+p collisions, the parameter T governs soft collisiona^l which is related to 
freeze out temperature for larger systems produced in relativistic heavy ion collisions 
after collective effects are separated out. The parameter n gives a good idea of initial 
production. A value of n close to 4 indicates point-like qq scattering (leading twist) 
while a large value of n is indicative of multiple scattering centers (higher twist 
effect).l^'iSl A successful description of particle spectra with the Tsallis distribution 
allow us to calculate the integrated yield which provide important information on 
the bulk properties of the soft particle production and also used to infer the degree 
of chemical equilibration when compared with the thermal modeh'^^ 

In the present work we use Tsallis distribution fit to draw systematics from the 
transverse momentum {pr) spectra of identified hadrons measured in p+p collisions 
at RHIC (Vs = 62.4 and 200 GeV) and at LHC energies (^/s = 0.9, 2.76 and 7.0 
TeV). We obtain the Tsallis parameters T and n for pions, kaons and protons and 
study them as a function of center of mass energy. 



2. Particle spectra and the Tsallis distribution 

The characteristics of transverse momentum spectra of particles produced in proton- 
proton or heavy ion collisions are keys to understand the particle production mech- 
anisms. Hagedorn,^ successfully described the shape of transverse mass rriT spectra 
of hadrons produced in heavy ion collisions in terms of two parameters by this form 
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Here A, pq and n are fit parameters. The parameter A is related to dN/dy which we 
will show little later. At low transverse momenta it assumes an exponential form and 
at large transverse momenta it becomes a power law which mimics "QCD inspired" 
quark interchange modePSl as follows: 



/, , mT\ " f-nmT\ 



exp , for pT -> (2) 

V Po 

~ ( — ) , for pT ^ oo. (3) 

The Hagcdorn function has been successfully used to describe the meson spectra 
in p+p collisions at RHIC energies with an additional parameter to describe heavy 
ion collisions as well.'ii' 

The Tsallis distributiorP'^l describes a thermal system in terms of two parameters 
T and q and is given by 

/ p\ -1/(9-1) 

G,iE) = cJl + iq-l)-) . (4) 



Here Cq is the normalization constant, E is the particle energy, T is the Tsallis 
temperature and q is the so-called nonextensivity parameter which measures the 
temperature fluctuations^] in the system as: q— 1 = Var{T)/ < T >^. The values of 
q lie between 1 < q < 4/3. For g — >■ 1, the distribution corresponds to an equilibrated 
system described by a pure exponential (Boltzmann-Gibbs) type distributionPSESl 

G,(i?) ^Cicxp(^-|^. (5) 

Using the relations E — niT — y/pj^+rn^ at mid-rapidity and l/(g— 1) = n, Eq. |4] 
takes the form 

which is same as Eq. [T] with po = nT. Larger values of q correspond to smaller 
values of n describe a system away from thermal equilibrium. In terms of QCD, 
smaller values of n imply dominant hard point-like scattering. Phenomenological 
studies suggest that, for quark-quark point scattering, n ^ A^^^aaA when multiple 
scattering centers are involved n grows larger and can go upto 20 for protons. 

In order to associate the Tsallis distribution with a probability distribution, 
which describes the invariant particle spectra defined over Q < E < oo, Eq. |6] 
must satisfy a normalization and energy conservation condition. Using the unit 
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normalization condition, we can determine the co-efiicient C„, the resulting formula 
used for fitting the hadron spectra used by PHENIX-21 and is given by 

d^iV _ I dN (n-l)(n-2) / nr + mr V" 

dp^ 2tt dy [nT + m{n ~ l)){nT + m) \ riT + m J 

In some of the papers,'^'^ another variant of Tsallis distribution is written as 

/ p\ -9/(9-1) 

G,iE)^C,i^l + {q-l)-j . (8) 

The corresponding normalized distribution in terms of n can be obtained by replac- 
ing n by (n — 1) in Eq. [T) 

The Tsalhs distribution used by some pubhcations of PHENIX,^^! STAEPIH] 
and CMy22l is of the form given by 

dp^ 2Tr dy {nC + m{n - 2)) (nC) \ nC J ^ '' 

This is same as Eq. [71 with the parameter C of Eq. [9] related to T by nC — > 
nT + TO. In summary, all the forms of Hagdorn and Tsallis distribution described 
above are essentially same with parameter of one form is related to the parameter 
of the other form by a definite relation. We will use the form given in Eq. [7] in rest 
of the analysis. 

3. Results and Discussions 

In the present work, the measured pt distribution of hadrons are taken from 
PHENIX {\y\ < 0.35),E3E2l|23l stAR (|j/| < 0.5, \y\ < 0.75)!^[2i] and CMS 
{\y\ < 1.0)221 experiments. The average yields of charged pions, charged kaons and 
protons measured at all RHIC {^/s = 62.4 and 200 GeV) and LHC energies {^/s = 
0.9, 2.76 and 7.0 TeV) are used in the analysis. The errors on the data are taken as 
quadratic sums of statistical and uncorrelated systematic errors wherever available. 

First, we fit all the measured hadron (pion, kaon and proton) spectra with 
Tsallis distribution (Eq. ^ keeping both the parameters T and n as free. Figure [1] 
shows the variation of parameters T and n as a function of ^/s along with their 
parameterizations by a function: a — 6(y^)~". We notice that the parameter T 
has a decreasing but weak dependence on ^/s and the parameter n has a definite 
decreasing trend with increasing y/s. In case of pions, both the parameters vary 
smoothly as a function of collision energy. The solid line in Fig. [Ija) is obtained by 
fitting n for all particles with a single parameterized form. The solid line in Fig. [1] 
(b) is the same for T for all particles. These curve mostly pass through the pion 
points. 

There is correlation between parameters T and n, if any of the two increases, the 
other also increases and thus we need to fix one of them by some method and study 
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the behaviour of the other as a function of particle type and -^s. The parameter 
T as defined in Eq. [7] can be assumed to be same for all particles at a particular 
energy. We can fix it to the value averaged over all particle types or simply use 
pion T for all the particles. Thus first we fit the measured pion spectra with Tsallis 
distribution by keeping parameters T and n as free and then fix T for kaon and 
proton spectra to get n. 
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Fig. 1. The variation of Tsallis parameters (a) n and (b) T as a function of center of mass energy 
{\/s) for pion, kaon and proton. The solid curve in each figure are the function a — b(^/s)~" . 

fitted over all particles. 
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Figure [2] shows the invariant yields of pions, kaons and protons as a function of 
PT for p+p coUisions at (a) ^/s = 62.4 GeVW^ and (b) ^/s = 200 GeVW^^ 
Figure |3] shows the invariant yields of pions, kaons and protons as a function oi pT 
for p+p coUisions at (a) y/s = 0.9 TeV,[22! (b) ^i = 2.76 TeV^ 

Figure |4] shows the invariant yields of pions, kaons and protons as a function of 
PT for p+p collisions at 7.0 TeVW^ 

The solid curves in all the figures are the fitted Tsallis Distribution of Eq. ([7|) . 
In case of pions both T and n are kept free during fits. For kaons and protons the 
values of T is fixed to the value obtained for pions. The values of T and n for pions 
and n for kaons and protons so obtained are given in Table [1] at different collision 
energies. It can be noted that that we get good x^ values for fits at all energies with 
such procedure. In Fig[5l upper plot (a) shows the variation of Tsallis parameter n 
and the lower plot (b) shows the variation of ratio riproton/npion and rikaon/npion as 
functions of ^/s. It is observed that the parameter n, for pions, kaons and protons, is 
monotonically decreasing with increasing -y/s, which could be understood in terms 
of changing production mechanism at different collision energies. 

The QCD cross sections scale as l/p^ in terms of the power n = 2 rta — 4 where 
Ua is the number of participating quarks.'^'^ If the dominant process in meson or 
baryon production is point like scattering qq — ^ qq (referred as leading twist) the 
number of participating quarks is 4 and hence n = 4. The power can go up due 
to contribution of higher twists processes explained in the following. The mesons 
scattering on quark qn — > qn gives n = 2x6 — 4 = 8. For proton production if the 
subprocess qq — >■ qqqq is dominant then n = 8. Proton scattering process on quark 
{qp — > qp) gives n — 12. The value of n for pions could go up to 16 for process 
pn -> pn and that for protons could be 20 corresponding for process pp — >■ pp. 

The decreasing value of parameter n, for pions, kaons and protons with collision 
energy in Table [T] indicates that the dominant production process is moving from 
higher twist to leading twist as one goes up in yi. Wc find that the separation of 
n among pions, kaons and protons decreases as we go from lower center of mass 
energy (at RHIC) to higher center of mass energy (at LHC). This indicates that at 
higher collision energies the production process is more like qq point scattering for 
all particles which gives lower and similar values of n for different particles. 

Lastly, we study all hadron spectra at all energies using fixed values of T. From 
the free fit results given in Fig[ljb) we choose two values T — 110 MeV and T=95 
MeV and fit all hadron spectra to obtain power n for pions, kaons and protons as 
a function of a/s. The results for T = 110 MeV and 95 MeV are given in Figs. El^a) 
and (b) respectively. We find that for both the values of T, the parameter n as a 
function of ^/s decreases which is noticed in all our analysis. Although the values 
of n at r = 110 MeV are greater than their respective values at T = 95 MeV the 
general conclusions about n remain the same. 
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4. Conclusion 



In the present work we use Tsallis fit to draw systematic trends from the transverse 
momentum spectra of identified hadrons measured in p+p collisions at RHIC and 
LHC energies. We also review various forms of Hagedorn and Tsallis distributions 
and show their equivalence. We obtain the Tsallis parameters T and n for pions, 
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Fig. 2. The invariant yields of pions, kaons and protons as a function of px at (a) y/s = 62.4 
GeV,EIlllSl and (b) y^ = 200 Ge\G3[20l[IIl for p+p collisions. The solid curves are the fitted Tsallis 
distributions. 
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Fig. 3. The invariant yields of pions, kaons and protons as a function of as a function of pT for 
p+p collision at (a) y/I = 0.9 TeVpl (b) ^ = 2.76 Te\^2I] and The solid curves are the fitted 
Tsallis distribution. 



kaons and protons and study them as a function of center of mass energy. Since T 
and n are correlated we fix the value of T for all particles at a particular energy to 
the pion temperature. In general, the Tsallis temperature T has a decreasing but 
weak dependence on center of mass energy which means there are less soft collision 
processes as we move up in energy. The power n determines if the particle is coming 
from point like qq scatterings (leading twist) or from multiple scattering centers 
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Fig. 4. The invariant yields of pions, kaons and protons as a function of as a function of pT for 
p+p collision at 7.0 TeV .1221 The solid curves are the fitted Tsallis distribution. 

involving many quarks (higher twists). The power n decreases with decreasing ^/s 
for all particles in all our analysis which indicates that the production process is 
moving from higher twist to leading twist as one goes up in yfs. The maximum 
values of n are found to be 16 for proton spectra 12 for pion spectra at the lowest 
RHIC energy. Another important observation is with increasing ^/s, separation of 
powers for protons and for pions narrows down hinting that the baryons and mesons 
are governed by same production process as one moves to the highest LHC energy. 
Separation of n for protons and pions is more at lower energy. This would imply 
that higher twist process prefer pion production over proton at lower energies. 
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Table 1. The values of Tsallis Parameters (T and n) for pions, kaons 
and protons. In case of pions both T and n are kept free during fits. 
For kaons and protons the values of T is fixed to the value obtained 
for pions. 





V^ (TcV) 


T (McV) 


n 


dy 


xVn^i/ 


For pion 


0.062 


117.6 ± 1.2 


11.87 ± 0.48 


0.812 ± 0.03 


0.19 


0.2 


112.6 ± 2.1 


9.30 ± 0.12 


9.204 ± 0.05 


4.52 


0.9 


104.5 ± 3.1 


7.16 ± 0.22 


1.917 ± 0.014 


1.18 


2.76 


97.5 ± 2.6 


6.08 ± 0.17 


2.44 ± 0.02 


1.04 


7.0 


94.4 ± 2.8 


5.55 ± 0.15 


3.083 ± 0.03 


0.88 


For Kaon 


0.062 


117.6 


10.66 ± 0.28 


0.072 ± 0.002 


0.17 


0.2 


112.6 


8.48 ± 0.13 


0.493 ± 0.008 


0.68 


0.9 


104.5 


6.53 ± 1.40 


0.238 ± 0.002 


0.07 


2.76 


97.5 


5.41 ± 0.08 


0.322 ± 0.000 


0.74 


7.0 


94.4 


4.70 ± 0.10 


0.418 ± 0.006 


0.33 


For proton 


0.062 


117.6 


15.51 ± 0.35 


0.007 ± 0.00 


0.35 


0.2 


112.6 


11.84 ± 0.14 


0.083 ± 0.003 


1.18 


0.9 


104.5 


8.95 ± 0.07 


0.105 ± 0.00 


1.06 


2.76 


97.5 


7.28 ± 0.06 


0.137 ± 0.001 


1.35 


7.0 


94.4 


6.07 ± 0.07 


0.177 ± 0.001 


1.05 



